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:^ ■ Abstract 

X 

y3 , We present a new method for the derivation of convolution identities for finite sums 

of products of Bernoulli numbers. Our approach is motivated by the role of these 
identities in quantum field theory and string theory. We first show that the Miki 
identity and the Faber-Pandharipande-Zagier (FPZ) identity are closely related, and 
give simple unified proofs which naturally yield a new Bernoulli number convolu- 
tion identity. We then generalize each of these three identities into new families of 
convolution identities depending on a continuous parameter. We rederive a cubic 
generalization of Miki's identity due to Gessel and obtain a new similar identity gen- 
eralizing the FPZ identity. The generalization of the method to the derivation of 
convolution identities of arbitrary order is outlined. We also describe an extension 
to identities which relate convolutions of Euler and Bernoulli numbers. 



1 Introduction: Convolution identities for Bernoulli numbers 

The Bernoulli numbers Bn are defined by the generating function pQ 

oo ^ 

b(x) = ^^ = yBn^ (1.1) 

n=0 

As is well known, the i?„ play an important role in combinatorics and number theory, and there 
exist many combinatorial identities involving these numbers j21 El- In the present paper, we 
will be concerned with the special case of convolution identities, which involve finite sums of 
products of Bernoulli numbers. The best known such convolution identity was found already by 
Euler (and independently by Ramanujan): 

Theorem 1.1 : (Euler; Ramanujan) For integer n > 2, 

Y,LljB2kB2n^2k = -i2n + l)B2n (1.2) 

The proof follows directly from the definition (jl.lj) . noting that the generating function b{x) 
satisfies 

bix)"^ = {l-x)b{x)-xb'{x) (1.3) 

Equivalently, H1.2|) follows by comparing the series expansions of either side of the trigonometric 
identity 

coth2x = 1 - (cothx)' (1.4) 

Many more such identities, involving folded sums of Bernoulli numbers, have been found since 
Euler 's work (see, e.g., jHElini). Most of them are similar to Euler 's identity (|1.2() in the sense 
that they involve -^ rather than B„ itself, as is the case already for the defining formula (|1.H) . 
Identities involving the Bn themselves, without the factorial denominator, are much rarer. One 
such identity was found by H. Miki in 1978. 



Theorem 1.2 : (Miki 7 ): For integer n > 2, 

E 02k02n-2k _ ST^ ^2fc-P2n-2fc / ^"- \ -"2n jt (-, e-\ 

^_^{2k){2n-2k)~ f-^^{2k){2n-2k)\2kl^ r,. ^n I -^j 

Here if,- denotes the ith harmonic number, 



H^ ^ iZ- (1-6) 



.=1^ 



As is well-known, the harmonic numbers can be alternatively expressed in terms of the digamma 
function 'ip{x) = r'(a;)/r(a;) and the Euler- Mascheroni constant 7: 



Hi = V'(i + l)+7 (1-7) 



Miki's original proof of the identity (|1.5j) identity is quite involved [Zj. A more elementary proof 
was given recently by I. M. Gessel [8 , using two different expressions for the Stirling numbers 
of the second kind. In Section 2.1 we present an even simpler proof, based on an appropriate 
generating function. 

In 1998 C. Faber and R. Pandharipande [5] found that certain conjectural relations between 
Hodge integrals in Gromov-Witten theory |in| lllj (see also J12j ) require the following identity 
to hold. 

Theorem 1.3 : (Faber and Pandharipande, with a proof by Zagier P|) For integer n > 2, 

E B2kB2n-2k _ 1 V^ B2kB2n-2k f'^'>^\ -^2n rr f-io) 

^^^ {2k){2n -2k) ~ nf^^ {2k) \2k) ^ n '""' ^ ^ ^ 

/I — 2"^^\ 
where B^ = ( ^„_i J^n (1.9) 

A proof of the Faber-Pandharipande-Zagier (FPZ) identity H1.8() was given by D. Zagier in an 
appendix to 9 . Note that the structure of the FPZ identity 1)1. 8() is similar to Miki's identity 
()1.5() . We later show (see Theorem 2.1) that this similarity is even more striking if Miki's identity 
is written in a slightly different form. 

Apart from pure mathematics, the Bernoulli numbers appear prominently in perturbative quan- 
tum field theory. This comes about at a very basic level [12]: perturbative loop calculations in 
quantum field theory generally involve traces of inverse powers of derivatives of functions defined 
on a circle. Since the spectrum of the ordinary derivative operator dp with periodic boundary 
conditions consists of the integer numbers, one has 

oo 

M^p^'^) ~ E fcl^ = C(2n) (1.10) 

fc=i 



But C,{2n) is related to the Bernoulli numbers through Euler's identity, 

B2n = (-l)"+^2-^C(2n) (1.11) 



The Bernoulli numbers also appear naturally in so-called "effective action" computations in 
quantum field theory, a field pioneered by W. Heisenberg and H. Euler, and V. Weisskopf |14j . 

In J15|I1HJ. the present authors found that Miki's identity arises naturally in a certain computa- 
tion in perturbative quantum field theory. Specifically, it arises in the course of the calculation 
of the two-loop effective Lagrangian for quantum electrodynamics in a constant background 
self-dual field. This calculation was done using two different integral representations. It turned 
out that both representations yield a result for the coefficients of the weak field expansion of this 
effective Lagrangian which involve a convolution of Bernoulli numbers, and that what is needed 
to show the equivalence of both results is precisely Miki's identity ()1.5)) . Since the two integral 
representations used are related by a simple coordinate transformation, this actually yields a 
new, and quite straightforward, proof of Miki's identity. This proof will be given in section [2Tl 
The simplicity of the approach presented here suggests a number of generalizations of these two 
identities, some of which are presented in sections |21 to |S1 Further generalizations are outlined 
in the conclusions. 



2 Simple generating function proofs of Miki and FPZ identities 

In this section we presented simple unified proofs of the Miki and FPZ identities based on 
generating functions. 



2.1 Generating function proof of Miki's identity 

Consider the generating function 



ipix) = ^p{x)—lnx + 



1 

2x 



(2.1) 



This function plays an important role in the quantum field theory computations in J15l I16j . 
From the asymptotic (large x) expansion of the digamma function PP it follows that 



ip{x) 



E B2k 1 
9i. T,2fc 



k=l 



2k X^ 



(2.2) 



Thus for the square of ip one finds 



°° 1 "~^ R R 
''J ~ 2^ T.2n 2^ ('9i-V9r7 - 71 



X 
n=2 k 



' {2k){2n - 2k) 



(2.3) 



Thus, we see that the [^(x)]^ is the generating function for the left-hand side of Miki's identity 
()1.5() . We prove Miki's identity by comparing ()2.3|) with the square of the following integral 
representation [see Eq 1.7.2 (25) in ITTj] of ^(x) : 



ds e ^^"^ ( coth s 



'tp{x) 

/o ^ s 

We break this comparison into three straightforward lemmas 
Lemma 2.1.1 : 

du 1-1 + 2 



(2.4) 



ip^x) 



dyye-^'^y I dn <| -1 + 2 ( cothy - i 











coth yu 

yu 



u 



( coth yu ) — I coth y I 

V v^J V yJ 



y{l-u) 
Proof : Squaring (|2.4|) we find 

Ujj(x)f = ds ds'e~^^(''+''')icothscoths'- (-coths' + -coths) H -} 

Jo Jo ^ ^s s' J ss' ) 

Now, using the trigonometric identity 

coth s coth s' = coth(s + s')( coth s + coth s'j — 1 , 

together with the symmetry s '^ s' ., and the transformation of variables ^ 

,s' 
y = s + s' 



s + s' 



(2.5) 



(2.6) 



(2.7) 



(2i 



^In quantum field theory terms, this change of variables corresponds to a change from a Feynman parameter 
integral (see, e.g., (18) 'l to a worldline parameter integral 19, . 



it is straightforward to show that [V'l^;)]^ can be rewritten as in Lemma 2.1.1. Note that the 
change of variables 1)2. 8|) introduces a Jacobian factor of y. D 

Our proof of Miki's identity follows by evaluating the asymptotic expansions of the integrals 
appearing on the right-hand side of (|2.5)) . 

Lemma 2.1.2 : 



2 rdyye-'^y [' dJcothy - -) fcothyn - J-) r. y J-y ^ B2A,_,, /2n\ 
Jo J^ \ ^ y)\ ^ yu) ^^x^-f^^{2k){2n-2k)\2kj 



Proof : The u integral is elementary: 



1 



1 / sinh y 
-m 



du ( coth yu 
V yu/ y V y 

After an integration by parts, the remaining y integral takes the form 

2x fdye-'^y In' {^^) 

The asymptotic expansion of the y integral is obtained using the Taylor expansion ^ 



(2.9) 



(2.10) 



KT)^^ 



°o o2fc-l 



2 D2k 2k 



' k{2k)\ 

which directly yields the result (|2.9jl after performing the y integration. D 
Lemma 2.1.3 : 



(2.11) 



(2.12) 



-2 r dyye-'^y [ du { —^ 
Jo Jo U(l- 



u ( coth yu ) — ( coth y 

yuj V y 



-^>^ ^2n 



n=l 



^2n. ^ 



(2.13) 



Proof : First, consider the u integral. We use the Taylor expansion of the coth function ^, 



' coth y = y^ 



',2k 



k=0 



nZK Tj 
f__D2k 



-y 



2k 



(2.14) 



for both cothyti and cothy. The u integral becomes elementary: 



du— 



ufcothyn- ^j - (cothy- i 



{l-u) 



E 



02n2 y 



2n,,2n-l /■! 



E 

n=X 



(2n)\ 

02n2 y 



du 



n=l ^-^'- -^0 

°° T^ o2?i„,2n-l 



,2n 



l-U 



(2n) 



-H, 



2n 



(2.15) 



Doing the y integral we obtain the result of Lemma 2.1.3. D 

Miki's identity (|1.5j) is then proved by comparing the results of Lemmas 2.1.1 - 2.1.3 with (|2.3)) . 

We conclude this section on Miki's identity by remarking that in the proof of Lemma 2.1.2, the 
partial integration in y leading to (|2.11|) is not essential. If, instead, one does the y integral 
directly using (|2.14|) and (|2.12|) . one arrives at a slightly different version of Miki's identity: 

Theorem. 2.1 : (Modified form of Miki's identity): For integer n > 2, 

E B2kB2n-2k _ j_ V^ B2kB2n~2k /^2n\ -B2n rr 

(2fc)(2n - 2k) ~ n^ i2k) y2k) ^T ^'^ 



^E^^ifrC:)-^^- (-«) 



where we have used H2n = -f^2n-i + ^, and Bq = 1. 

Comment 2.1.1 : This last form 1)2. 16j) of Miki's identity brings out most clearly the similarity 
to the FPZ identity (fTH)) . 

2.2 Generating function proof of the FPZ identity 

To prove the FPZ identity we use, instead of tpix), the generating function "iplx) defined by 

^5(3;) = 'iIj{x + -) -Inx (2.17) 

The large x expansion of ip{x) is 

fc=l 

where i?2fc was defined in (|1.9|) . The expansion 1)2. 18() follows from the corresponding expansion 
p.2j) for ip{x), using the "doubling" identity ^ for the ip function, 

^(2x) = -V'(x) + -V'(x + -) + ln2 (2.19) 

Thus, the square of 'tp{x) is the generating function for the left-hand side of the FPZ identity 

The generating function for the right-hand side of the FPZ identity is obtained by squaring the 
following integral representation for ip{x): 

i(.) = -l «»,--(_--) (2.21) 



Lemma 2.2.1 : 



dyye ^^^ du < ( cothyn 

Jo I smh y 



1 



2/(1 -u) 



u 



1 



sinh yu 



1 

yu 



1 



1 

yu 
1 



sinh y y 



Proof : Squaring the integral representation (|2.21|) . using the trigonometric identity 

1 _ coth(s) + coth(s') 

sinh(s) sinh(s') sinh(s + s') 

and the symmetry under s -^-^ s', it follows that 

ds'e-2( 

1 1 



(2.22) 



(2.23) 



m^)Y 



2 / ds 
'o ^0 

1 

+ - 



'^'i^^ -fcoth(.')-^ 

\ sinh(s + s')\ ^ ' s' 



sVsinh(s + s') s + s' sinh(s') s' 
Applying the transformation of variables (|2.8|) , we obtain (|2.22)) . D 



(2.24) 



Our proof of the FPZ identity now follows by evaluating the asymptotic expansions of the 
integrals appearing on the right-hand side of 1)2. 22() . 



Lemma 2.2.2 



dyye ^"^^ / du 
Jo 



sinhy 



coth yu 



1 

yu 



t^i±^^('") p-^^) 



2k 



Proof : The proof is almost identical to the proof of Lemma 2.1.2, but in doing the y integral 
we use an asymptotic expansion of 1/ sinh y rather than coth y. This has the effect of replacing 
one of the Bernoulli number factors i?2n-2fc by B2n-2kj and also of changing the upper limit of 
the k summation from in — 1) to n. D 



Lemma 2.2.3 



-2 / dyye-^'-'y / du 
'0 Jo 

oo 



1 



y(l -u) 



u 



1 



sinh yu 



1 

yu 



sinhy 



E 

n=l 



1 B' 



2n 



^2n ^ 



H, 



2n~l 



(2.26) 



Proof : The proof is almost identical to the proof of Lemma 2.1.3, except we use the asymptotic 
expansion of 1/ sinhy rather than (cothy — -), which has the effect of replacing i?2n by B2n- D 

The FPZ identity (|1.8|) is then proved by comparing the results of Lemmas 2.2.1 - 2.2.3 with 
()2.2U|) . Note that the n = 1 terms on the right-hand sides of ()2.25|) and (|2.26|) cancel, permitting 
the comparison with ()2.2fl|) . 

3 A new convolution identity 



The similarity between the proofs and forms of the Miki and FPZ identities immediately suggests 
a new identity, in which on the left-hand side the B2n and B2n are mixed. As is clear from Section 



2, such an identity could be derived using the generating function ip{x) ^(x), and comparing its 
summation and integral representations. However, there is another, even simpler, way to derive 
this mixed identity. Note that the two generating functions ip^x) and ^(x) are related via the 
■0 function doubling identity (|2.19j) as: 

ip{x)+ij{x) = 2'ijj{2x) (3.1) 

Thus, it follows that 

r ~\ 2 r ~\ 2 

2'ijj{x)'4>{x) = 4: tl){2x) - ^{x) -[i>{x)Y (3.2) 

A new Bernoulli convolution identity emerges by using the asymptotic expansions (|2.2)) and 
()2.18|) for the left-hand side of ()3.2() , and the asymptotic expansions of the squares of the integral 
representations (|2.4j) and (|2.2H) of ip{x) and '^(x), respectively, on the right-hand side. All 
necessary results for the squares of the relevant integral representations are contained in Lemmas 
2.1.1 - 2.1.3, and 2.2.1 - 2.2. 3. To express the result in a symmetrical form, we use the modified 
form of Miki's identity in ()2.16|) . 

Theorem 3.1 : For integer n > 2: 

BokBon-oh 1 ^-^ B').kB2'n-').k f2n\ /I — 2 ~\ IB 



2k02n-2k _ i_ V^ -P2fc-P2n- 2fc ( ^n \ ( I — Z, \ I 02n 

' {2k){2n -2k) ~ n^ {2k) 



E J^2k^2n-2k _ ^ ST^ -"2fc-"2n-2fc I ^n \ i L — ^ \ i U2n tt (n o\ 



Proof : The generating function for the left-hand side is given by half the left-hand side of (|3.2)) . 
The right-hand side is obtained by using Miki's identity in the form (|2.16() for the squares of ip, 
and the FPZ identity (|1.8j) for the square of '(/'■ Simple algebra then leads to the form in (|3.3|) . 



4 Three infinite families of convolution identities 

The use of the generating functions ip{x) and i/)(x) to prove the Miki and FPZ identities, as well 
as the new "crossed" identity ()3.3|1 in Theorem 3.1, immediately leads to natural generalizations 
of each type of identity. 

4.1 Generalization of Miki's Identity 

To derive a generalization of Miki's identity, consider the p^^ derivative of the generating function 
ip{x). This has the large x asymptotic expansion: 

^ ^^ ^ ' ^ (2n)r(2n) x^^+p ^ ' 

This function also has the following integral representation: 

^(P)(x) = -(-2)P / dse-^'^'sP (^coths - -j (4.2) 

We can use this integral representation to extend the definition of ^p^^' {x) to non-integer values 
of p; for this extrapolation the expansion (|4.H) continues to hold, as can be seen by using (|2.14j) 



under the integral in (|4.2|) . Thus, in the following let p denote an arbitrary non- negative number. 
We can derive new identities, for any such p, by squaring these two representations of ijj^P'{x), 
and then comparing, just as was done (for p = 0) to prove Miki's identity. The proof proceeds 
in a very similar manner. 



Lemma 4.1.1 : 

n-l 



J.fo)^^^f ■ V 1 V B,,B,^^,, T{2k + p)n2n-2k+p) 

J ^a;2"+2p^ (2A;)(2n-2A:) r(2A:)r(2n - 2A:) ^'' 

Proof: follows from (jUTJ. D 
Lemma 4.1.2 : 

V'(P)(x)] ^ = 22P r dyy^P+^ e-^'^y f duuP{l - nf |-1 + 2 fcothy - -^ (cothyu - —^ 



2 



y{l-u) 



coth yu ) — I coth y I 

yuj V yJ 



(4.4) 



Proof : Square the integral representation (|4.2j) . change variables from s and s' to y and u, as 
in ()2.8|) . and regroup terms as in the proof of Lemma 2.1.1. Note that the argument about 
symmetrizing with respect to s and s' still holds because the extra factors in the integrand 
appear as {s s'Y = y^^u^(l — u)^. D 

Now consider each of the three terms appearing on the RHS of ()4.4() . 

Lemma 4.1.3 : 

- 2^^ j^ dy /^+^ e-'^y [duu^il- uf ~ -^-^^^ (4.5) 

Proof : immediate. D 



Lemma 4.1.4 : 



22P+1 r dyy^P+^ e-^'^y f duuP{l - uf (cothy- -^ (cothyu - —^ 

9rr I ^)Y^ ^ V B2kB2n-2k T{2k + p)T{2n + 2p) 

{P+ ) Z^ ^2n+2p Z^ (2kM(2n - 2k)\ r(2n^2k^l) ^ ' 



n— 1 

^_2 x""-^"^ ^ (2A;)!(2n - 2k)l T{2p + 2k + 1) 
Proof : First, consider the u integral: 

"1 / 1 \ °° n o2n„,2n-l /■! 



fl / 1 \ R r)2n„2n~l /■! 

,22V""^r(p + i)r(p + 2 

(2n)! r(2p + 2n + l 



n)P 



n=l 
°° " o2n,,2n-l 



^ i?2n2""j/^"-^ T{p + l)r(p + 2n) 



n=l 



Now doing the y integral we obtain: 



^ (2n)! r(2p + 2ra + l) Jq \ V 

g2p+i V ^2n22" T{p + l)r(p + 2n) ^ B2fc2=^fe r(2p + 2n + 2fc) 

^ (2n)! r(2p + 2n + l) ^ (2fc)! (2x)^P+2n+2k ^ ■^) 



from which Lemma 4.1.4 fohows. D 
Lemma 4.1.5 : 

1 



1 



Jo Jo U(l-^i) 

~ 2^:;2^T2? j^^i '^Pip + k,p + i) 



u I coth yu 



I coth y 

V y 



(2n)! 



(4.9) 



n=2 " ^ ' k=l 

Proof : First, consider the u integrah 

•1 



duvF{l — u) 



p-i 







E 

n=l 

oo 

E 



^2n^ y 



u I coth 

2n„,2n-l />1 



yu 



yu 



I coth I 



\p-l (^,2n 



(2n)! JO 

2n„,2?i-l 2n 



^2n^ y 



(2n)! 



J^/3(p + /c,p+l) 



(4.10) 



where /3(p, q) is the Euler beta function. Doing the y integral we obtain the result of Lemma 
4.1.5. D 

We are now ready to state the generalization of Miki's identity: 

Theorem. 4.1 : For any p >0, and for integer n > 2: 

^' ^2fc-B2n-2fc r(2fc+p)r(2n-2fc+p) ^ +nV ^^^^^^-^^ T{2k+p)T{2n + 2p) 



^{2k){2n-2k) T {2k)T {2n - 2k) '-^ ' ^ {2k)\{2n - 2k)] T{2p + 2k + l) 



fc=i 



Proof : Follows by comparing the result of Lemma 4.1.1 with those of Lemmas 4.1.2 - 4.1.5. D 

Comment 4.1.1 : When p = we recover from Theorem 4.1 Miki's identity in the form of 
Theorem 2.1. 

Comment 4.1.2 : When p = 1 we obtain from Theorem 4.1 a convolution identity which just 
involves the Bernoulli numbers themselves on the left-hand side: for n > 2, 



Y^ B2kB2n^2k - —— Y^ B2kB2n-2k ( 2fc , 2 ) 
fc=l fc=l \ "T / 



2n + 2\ 

2n-2k\ „, , o + 2ni?2ji 



(4.12) 



4.2 Generalization of the FPZ Identity 

To derive a generalization of the FPZ identity, consider the p^^ derivative of the generating 
function 'il>{x). This has the large x asymptotic expansion: 



,{;(P)ix)^i-l)P+'Yl 



B2nTi2n + p) 1 



^ (2n)r(2n) a;2"+P 



(4.13) 



10 



This function also has the following integral representation: 



^(p)(x) = -(-2y 



dse-^'^'sP 



1 



sinhs 



(4.14) 



As in the Miki case, we can use (|4.14|) to define ^'^^'(x) for non-integer p. We can then derive 
new identities, for any positive p, by squaring these two representations of 'ip^^'{x), and then 
comparing, just as was done (for p = 0) to prove the FPZ identity. The proof proceeds in a very 
similar manner. 



Lemma 4.2.1 



n-l 



il:^P\x) 



2-^ ni^2n+2p 2-^ 



B2kB2n^2k r{2k + p)T{2n - 2k + p) 



ra=2 



2n+2p z_> (2k){2n - 2k) T{2k)T{2n - 2k) 



(4.15) 



Proof: follows from (|I?T^ . D 
Lemma 4.2.2 : 



ij^-P^x) 



22p+l 



dy y2p+i e'2^^ / du uP(l - u)p | ^— f coth yu - — 

[ sinh y \ yu 



1 



y[l-u) 



e 
1 



1 



1 



1 



(4.16) 



^sinhyn yuj \sinhy y^ 

Proof : Square the integral representation (|4.14)) . change variables from s and s' to y and n, and 
regroup terms as in the proof of Lemma 2.2.1. D 

Now consider each of the two terms appearing on the RHS of H4.16|l . 

Lemma 4.2.3 : 

22P+1 r dyy^P-^^e-^'^y I duu^il - uY [-^ (cothyu - - 
Jo Jo [smhy \ yi 



2rf , 1^^p 1 V B2kB2n-2k T {2k + p)T {2n + 2p) 
^P+ > 2^ ^2n+2p 2^ ^2k)\{2n - 2k)l r{2p + 2k + l) 



(4.17) 



n=l 



Proof : The proof is almost identical to the proof of Lemma 4.1.4, but in doing the y integral we 
use an asymptotic expansion of l/sinhy rather than cothy — -. This has the effect of replacing 

one of the Bernoulli number factors i?2n-2fc by B2n-2k^ and also of changing the upper limit of 
the k summation from (n — 1) to n. D 

Lemma 4.2.4 : 



Jo Jo ' {y{l-u) 



1 



u 



1 



1 



sinh yu yu 



1 



1 



sinh y y 



2n 



^^ 1 B2nT{2n + 2p) ^^^ 
^2:^Z2^^T2^ 7^^ 2^P{p + k,p + l) 

n=2 ^ '' k=\ 



(4.18) 



Proof : The proof is almost identical to the proof of Lemma 4.1.5, except we use the asymptotic 
expansion of 1/sinhy rather than cothy. D 

We are now ready to state the generalization of the FPZ identity: 



11 



Theorem 4.2 : For any p > 0, and for integer n > 2: 

y.^ B2kB2n-2k T {2k + p)T {2n - 2k + p) ^ y. B2kB2n-2k T {2k + p)T {2n + 2p) 

^{2k){2n-2k) T{2k)T{2n-2k) ^^ ' ^ {2k)\{2n - 2k)\ T{2p + 2k + l) 

+2 ^-""g"+^^' 'i:'/;(p+>..p+i) (4,19) 



fc=i 



Proof : Follows by comparing the result of Lemma 4.2.1 with those of Lemmas 4.2.2 - 4.2.4. D 

Comment 4.2.1 : When p = we recover from Theorem 4.2 the FPZ identity (|1.8)) . 

Comment 4.2.2 : When p = 1 we obtain from Theorem 4.2 a convolution identity which just 
involves the Bernoulli numbers themselves on the left-hand side: for n > 1, 



^ - - 1 V^ D D /2n + 2\ - 

/ ^ -D2fc-D2n-2fc = — — j- 2_^ B2kB2n-2k 1 r,r, , o ) + 2ni?2n 

k=i "■ "*" fc=i V + / 



(4.20) 



4.3 Generalization of Theorem 3.1 

To generalize Theorem 3.1, we differentiate p times the relation (|ci.lj) connecting the two gener- 
ating functions 'ip{x) and ip{x). This leads to 



^(p)(a;) +^(p)(x) = 2P+i^(P)(2x) 

This relation also holds true for non-integer positive p, as can be easily seen using the 
representations ()4.2() . 1)4.14(1 for tp^P' and tjj^P' and the trigonometric identity 



sinh(s) 



coth(; 



coth(s) + 

Squaring the relation (|4.21l) we obtain 
Lemma 4.3.1 : 

2ip'^P\x)^p'-P\x) = 22P+2 U(p)(2a;)j^ - \i;'^P\x)\^ - \^^p\x) 

{p >0). This brings us to 

Theorem 4.3 : For any p >0, and for integer n > 2: 

y.^ B2kB2n-2k r(2fc + p)r(2n -2k+p) _ 
^ {2k){2n - 2k) T{2k)T{2n - 2k) 

5,i.B9„_,i. /l-22^-i\r(2A: + p)r(2n + 2p) 



2r(p + i)^ 



^ {2k)\{2n - 2k)\ 



I \ 2^"^^ 



r(2p + 2A; + l) 



2n-l 



^^liy3?^E«.^'-^') 



(2n)!22 



fc=i 



(4.21) 
integral 

(4.22) 



(4.23) 



(4.24) 
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Proof : The proof follows by taking the product of the expansions of the two functions on the 
left-hand side of Lemma 4.3.1, and comparing with the expansions of the squares of the integral 
representations of the three terms appearing on the right-hand side of Lemma 4.3.1, using the 
results of Theorem 4.1 and Theorem 4.2. D 

Comment 4.3.1 : When p = we recover the identity in Theorem 3. 1. 

Comment 4.3.2 : When p = 1 we obtain from Theorem 4.3 a, convolution identity which just 
involves the Bernoulli numbers themselves on the left-hand side: for n > 1, 

"-1 _ 1 " /I - 2^^^^ /2n + 2\ B 

Y: B2kB2n-2k = — Y E B2kB2n-2k (^l^^T^) Ufc + 2 + ^^^ " ^^ 2^ ^^'^^^ 

k=l fc=l \ "T / 

Comment 4.3.3 : Note that in all the above the positiveness condition on p was used only to 
avoid singularities. Theorems 4.1, 4.2 and 4.3 actually hold true also for negative p as long as 
none of the T - factors on either side becomes singular. 



5 Higher order convolution identities 

In the recent |S], I.M. Gessel shows the existence of an infinite tower of convolution identities 
involving multiple products of Bernoulli numbers, of which Miki's identity (|1.5|) is just the lowest 
order one. He also explicitly obtains the next element of this series, a triple product identity: 

Theorem 5.1 : (8', eq.(4)): For integer n > 3, 



E B2kB2lB2m _ sr^ B2kB2lB2m f 2n \ ,oTT Sr^ ('^^\ B2kB2n-2k 

^^_ (2A;)(2/)(2m) " ^^^^J2k){2l){2m)\2k,2l,2m)^ '" ^ W (2A:)(2n - 2^ 



k-\-l-\-m,^n 



Here 



k,l,m>l 



. „jj B2n / 2 3 5\ B2n-2 ,^^. 

+bHon 2 [n n-\ — ^ 5.1 

' 2n \ 2 4/ (2n-2) ^ ' 



l<«<i<2ra -^ ^ i=l ^ 

In our present approach, it is clear how to construct generalizations of the Miki and FPZ 
identities involving N - fold products of Bernoulli numbers: 

1. Take the A^th power of iIj{x) (resp. '4^{x)). The expansions (|2.2j) (resp. (|2.18j) ) generate 
the A^ — 1 - fold convolution on the left-hand side of the identity, 

Ei E nlr (^-3) 

■^i — 1 * 

(with Bi replaced by Bi in the FPZ case). 
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2. Use ((231) (resp. ^23l) to rewrite 

(V;(x))^ = / dsi ds2--- ds7ve-2^(^i+^2+...+siv)-QJ^coth(si)-- 

(5.4) 
(with coth replaced by ^^ in the FPZ case). 

3. Use trigonometric identities together with the symmetry of the integrand under permu- 
tations of the variables {si, . . . , sn} to rewrite the integrand in such a way that only the 
combinations si + S2 + • • • + sn, S2 + S3 + . . . + sn, • • • , sn-i + sn, sn appear as arguments 
of trigonometric functions (it is easy to see that this is always possible). 

4. Perform the transformation of variables from {si, . . . , sn} to {y, ui,U2, ■ ■ ■ , un-i} where 

y = S1 + S2 + ■■■ + SN 

SM+l + SM+2 + ■■■ + SN nj 1 AT ^ ^k k\ 

UM = ■ \ ■ , M = l,...,iV-l (5.5) 

Sl + S2 + ■ ■ ■ + SN 

The Jacobi factor of this transformation is y ■ 

5. Use the Taylor expansions ()2.14() . H2.12|) to do all integrals. 

Let us carry this through explicitly for the case A^ = 3. In the Miki case, after step 3 one finds 

Lemma 5.1: 

00 roo roo 



(,P(x)f = / dsi dS2 / dS3e'2x(si+S2+S3) 

Jo Jo Jo 



X<!6Ci23C23C3 + - 



^ ■S23 ^ 

<-^123 (-^23 

S123 



C3 H C123 

S2 



6 
+— 

S2 



— {Ci23 ^23) (C'123 ^^3) 

Si Si23 S12 S123 



(^23 L/3 

S23 

3C3 — 2C123 



S123 



(5.6) 



Here for compactness we have introduced the abbreviations 



S23 = 'S2 + S3, S123 = Si + S2 + S3, C(.) = COth(s(.)) . (5.7) 

S(.) 

Moreover, we have already combined terms in a way which facilitates the evaluation of the 
integrals (in particular, it avoids the appearance of spurious singularities). 

After the transformation 1)5. 5(1 the integrals can be done in a way which is completely analogous 
to the N = 2 case treated in section |21 The result is a slight modified form of the identity (|5.1|) , 

Theorem 5.2 : (Modified form of Gessel's identity): For integer n > 3, 

'2n\ B2kB2n~2k 



V^ B2kB2lB2m _ _3_ >r^ B2kB2lB2m f 2n \,'^tt S^{'^^\ 

{^_ {2k){2l){2m) - "^.^{^^ {2k){2l) {2k,2l,2m) n '""^^W) 



k-\-l-\-'m^n 



k.l,m>l 



m 



+QH2n 2^^ - n'--n + -) r'''-\ (5.8) 

' 2n V 2 4;(2n-2) ^ ' 
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B2n / 2 3 5\ B2n-2 



Comment 5.1 : The somewhat different form of the right hand side compared to (|5.1|) is due to 
the same type of ambiguity (regarding integrations-by-parts) which was mentioned aheady at 
the end of section 2.1. 

The FPZ case is again similar, though shghtly more comphcated. After step three one finds that 
the integrand can be written in the fohowing way, 



Lemma 5.2: 






{i;{x)f -- 


/•OO /"OO /"OO 

= / dsi / ds2 / 
Jo Jo Jo 




l Si 




6 

+— 

S2 


\^((^ S23 

(,'^123 
Lsi Si23 



ds3e-2^('i+"2+^3) 



e '523 e 

'-'123 '-'23 

Sl23 



C3 H 'S'123 

S2 



-'S'23) ('S'123 'S'3) 

S12 S123 

6 6 

H [C23 — 'S'23]C'3 H [C23 — S23 — C3 + 53] 

S123 ■S2S123 



<-'23 ^3 

S23 
A 
25123 



S123 



(5.9) 



Here we have used one more abbreviation. 



^ sinh(s(.)) S(.) 



Theorem 5.3 : (Cubic generahzation of the FPZ identity): For integer re > 3 

n-l 



2re\ B2kB2n-2k 



E B2kB2lB2m _ _3_ ■^ B2kB2lB2m f 2re \ ,^TT \^ f'^^\ 
^^_^ (2fc)(2/)(2rre) " 2^^^^^^ {2k){2l) \2k,2l,2m) ^ n '" ^^ W {2k) 

k,l,m'>l 

+ ^ E (2 J ^(^2n-2fc - B2n-2k) + ^i/2n-l(i?2n " B 



fc+i + m— n 
fc,Z,m>l k,l,m'>l 

n-l 



2n — 02n 

fc=l 

+6^2n,2^ - ^^^^2n-2 (5.11) 

ire 4 

Note again the similarity with the Miki case, 1)5. 8(1 . 

Comment 5.2 : It would be straightforward to extend both H5.8() and (|5.11j) to continuous families 
of identities along the lines of section 4. 

6 Conclusions 

The method presented here clearly allows one to derive, with relatively little effort, convolution 
identities for Bernoulli numbers of the quadratic type as well as higher order ones. Clearly 
we have not been able here to explore all its ramifications. For example, it should be possible 
to derive "mixed" identities such as (|3.H|) also at the cubic or higher level. Another possible 
direction is to use other generating functions to generate related identities involving the Euler 
numbers. The simplest such case comes from considering the generating function 



r-oo 

-2xs , 



g{x) = / dse ^^*sechs 
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oo 



Then it follows that 



from which one finds 



n=0 ^ ^ 



r / m9 /"°° , 9„,, Incoshu 

y(.)f = 2y^ *»-""-sr7^ . (6.2) 



EP 17 2 ^-^ B2kB2n-2k f„2k A r,2fc-l /i r,2n-2fc-A / ^""A /« Q^ 

i?2.-2i?2n-2. = - ^ ^ (2 -lj2 (1-2 j(^2A:J (^-^^ 

fc=l A;=l ^ ^ 

(n > 1). Clearly one can generate many other such identities relating convolutions of Euler 
numbers to convolutions of Bernoulli numbers. 

We conclude by emphasizing again that the types of generating functions and identities discussed 
here show up naturally in perturbative quantum field theory computations at the second-order 
(or "two loop" ) level ^1 El HU] . We expect related multiple convolution identities of higher 
order to play a similar role for higher-loop contributions to the eff'ective Lagrangian in quan- 
tum electrodynamics beyond the two-loop level. The higher order FPZ type identities might 
correspond to new relations between Hodge integrals and thus be of relevance for topological 
quantum field theory and string theory. 

Ackno'wledgements: We are very grateful to Richard Stanley for correspondence, and to 
Albrecht Klemm for discussions. C.S. thanks the Institut des Hautes Etudes Scientifiques, Bures- 
sur-Yvette, and the Albert-Einstein-Institut, Potsdam, for hospitality. G.D. thanks the US 
Department of Energy for support under grant DE-FG02-92ER40716, and thanks the Rockefeller 
Foundation for a Bellagio Residency Award. We also acknowledge the support of NSF and 
CONACyT for a US-Mexico collaborative research grant NSF-INT-0122615. 



References 

[1] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, Dover, New York, 
1972. 

[2] L. Comtet, Advanced Combinatorics, Dordrecht, Boston, 1974. 

[3] T. M. Apostol, Introduction to Analytic Number Theory, Springer, New York, 1998. 

[4] K. Dilcher, "Sums of Products of Bernoulli numbers", J. Number Theory 60 (1996) 23. 

[5] I.-C. Huang and S.-Y. Huang, "Bernoulli numbers and polynomials via residues" , J. Number 
Theory 76 (1999) 178. 

[6] A. Basu and T. M. Apostol, "A new method for investigating Euler sums", Ramanujan J. 
4 (2000) 397. 

[7] H. Miki, "A relation between Bernoulh numbers", J. Number Theory 10 (1978) 297. 

[8] I.M. Gessel, "On Miki's identity for Bernoulli numbers", J. Number Theory, to appear. 



16 



[9] C. Faber and R. Pandharipande, "Hodge integrals and Gromov-Witten theory", Invent. 
Math. 139 (2000) 137, math.A G/9810173 , 

[10] T. Eguchi, K. Hori, and C.-S. Xiong, "Quantum cohomology and Virasoro algebra", Phys. 
Lett. B 402 (1997) 71. 

[11] E. Getzler and R. Pandharipande, " Virasoro constraints and the Chern classes of the Hodge 
bundle", Nucl. Phys. B 530 (1998) 701, math.AG/9805114 

[12] C. Liu, K. Liu, and J. Zhou, "Marino- Vafa formula and Hodge integral identities", 
math.AG/0308015 

[13] M.G. Schmidt and C. Schubert, "On the calculation of effective actions by string methods", 
Phys. Lett. B318 (1993) 438, arXiv:hep-th/9309055 . 

[14] For a review, see : G. V. Dunne, "Heisenberg-Euler Effective Lagrangians: Basics and Ex- 
tensions", arXi v:hep-th/0406216 , to appear in Ian Kogan Memorial Volume, M. Shifman 
et al (Eds.), (World Scientific, Singapore). 

[15] G.V. Dunne and C. Schubert, "Two-loop self-dual Euler-Heisenberg Lagrangians 
(I): real part and h elicity amphtudes", J. High Energy Phys. 0208, 053 (2002) 
|arXiv:hep-th/02050"04l . 

[16] G.V. Dunne and C. Schubert, "Two-loop self-dual Euler-Heisenberg Lagrangians 
(II): imaginary part and Borel analysis", J. High Energy Phys. 0206, 042 (2002) 
,arXiv:hep-th/0205005, . 

[17] A. Erdelyi (ed.). Higher Transcendental Functions, Vol. I, Kreiger, Florida, 1981. 

[18] C. Itzykson and J. Zuber, Quantum field theory, McGraw-Hill, New York, 1985. 

[19] M.G. Schmidt and C. Schubert, "Worldline Green functions for multiloop diagrams", Phys. 
Lett. B 331 (1994) 69, i arXiv:hep-th/9403158| . 

[20] G. V. Dunne, "Two-loop diagrammatics in a self-dual background," J. High Energy Phys. 
0402, 013 (2004) |arXiv:hep-th/0311167| . 



17 



